Mean-field quantum phase transition in graphene and in general gapless systems 
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We study the quantum critical properties of antiferromagnetism in graphene at T — within 
mean-field (MF) theory. The resulting exponents differ from the conventional MF exponents, de- 
scribing finite temperature transitions. Motivated by this, we have developed the MF theory of 
general gapless phases with density of states p(e) ~ |e| r , r > — 1, with the interaction as con- 
trol parameter. For r > 2, the conventional MF exponents a la Landau are recovered, while for 
— 1 < r < 2, the exponents vary significantly with r. The critical interaction is finite for r > 0, 
therefore no weak-coupling solution exists in this range. This generalizes the results on quantum 
criticality of the gapless Kondo systems to bulk correlated phases. 
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Graphene, a single sheet of carbon atoms, arranged in 
a hexagonal lattice provided us with lots of excitement 
recently [1(. Its novel electronic properties arise from the 
linear energy-momentum dispersion of electrons at low 
energies, forming a Dirac cone, resulting in a linearly 
vanishing density of states (DOS) around half filling, sim- 
ilarly to d-wave superconductors. It explains the semi- 
metallic response of graphene, placing it somewhere be- 
tween good metals with a constant DOS and insulators 
with a clean gap. Although most of its electronic and 
transport properties can be understood by analyzing the 
Dirac equation under appropriate conditions, the effect 
of electron-electron interaction on its Dirac quasiparticlcs 
remains a bit mysterious [2|. No conclusive evidence for 
any kind of phase transition of either pristine or doped 
graphene has been reported to date, and in spite of the 
increasing amount of theoretical work[3rjl3l|. no consen- 
sus is reached on the role of interactions. 

In addition to graphene, other systems with semi- 
metallic behaviour are also available. In the graphene 
family, trilayer graphene with appropriate stacking ex- 
hibits a diverging DOS[14j, |l5| as ~ jej" 1 / 3 . Around 



a two dimensional semi-Dirac point, such as the one in 
V02/Ti02 nano-heterostructures j 1 6J . the DOS varies as 
~ \f\e\. Disorder belonging to specific classes can mod- 
ify the exponent of the DOS of Dirac electrons to almost 
arbitrary power law 17j as ~ |e| r with < r < 1. These 



can also be generalized to three dimensions [18|, provid- 
ing us with a whole zoo of systems exhibiting a power-law 
DOS around the Fermi energy. 

Systems with gapless DOS (p(e) ~ \e\ r with r > — 1) 
are expected to modify the behaviour of correlated phases 
and alter their critical behaviour. While quantum impu- 
rity models (i.e. Kondo model) in metallic hosts with 
a power law DOS have been extensively investigatedllfj 
following the pioneering work of Withoff and Fradkm[20J , 
much less is known about the interplay of bulk corre- 
lations (e.g. Hubbard model) and power law DOS ex- 



cept for isolated efforts on specific models. On general 
grounds, the power-law vanishing DOS implies restricted 
phase space for quasiparticlcs, which should weaken the 
effect of interactions, analogously to low-dimensional sys- 
tems (i.e. phase transitions are less likely). 

This motivates us to study the critical properties of 
bulk correlated systems with a power law DOS within the 
mean field (MF) theory. Our interest is twofold: on the 
one hand, we want to understand how and under which 
conditions correlations can alter dramatically the physi- 
cal properties of graphene (e.g. by opening a gap in the 
spectrum) at T = 0; on the other hand, the quantum crit- 
ical properties of general bulk gapless phases have never 
been systematically studied (except for the Kondo prob- 
lem, which is local) and our work intends to fill this gap. 

We start with the Hamiltonian of graphene in tight- 
binding approximation, written as 
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where t is the hopping integral and /(k) = 1 + e" 
e ika 2 w }t n the lattice vectors ai and &2, a^ and &ker 
arc the annihilation operators of the Bloch states corre- 
sponding to the A and B sublattices, respectively [1(. The 
electron spectrum vanishes linearly at the Dirac points as 
£o(k) = zkhvpk at the corners of the Brillouin zone (K, 
K'), vf is the Fermi velocity. The DOS is linear in energy 
p(e) — 2\e\/(iTh 2 v F ). Since it differs significantly from 
that in normal metals (i.e. p(e) ~const), it is natural to 
ask how the interplay of interactions and the above linear 
DOS can modify the critical dynamics, which serves as 
the basic motivation for our study. For the sake of sim- 
plicity, we focus on antiferromagnetism within MF ap- 
proximation in the presence of a simple repulsive on-site 
Hubbard interaction, so that s = (%2rs(R))/N is the 
expectation value of the staggered magnetization where 
■s(R-) = 2<r f J (a^ cr a RCT -6^ cr 6R cr )/4 and TV the number of 



unit cells. The MF theory for the Hubbard model gives 
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which describes quasi-particles with the dispersion rela- 
tion e(k) = ±y / sl(h) + A 2 , A = Us is the energy gap 
in the spectrum. We introduce a cut-off energy W, pre- 
serving the total number of states as W — yy/3n\t\. 

The ground state energy of the half-filled system (pris- 
tine graphene) at T = reads as 



E(s) U(l + As 2 ) 4C/ 3 
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subject to minimization with respect to s. As opposed 
to standard MF results, there is no non-trivial weak cou- 
pling solution, the normal phase with s = extends until 
U < W. For bigger values of U, we enter into the antifcr- 
romagnetic phase with \s\ = (1 — W 2 /U 2 )/2, the critical 
value of U being U c = W « 2.33|i|. The accurate MF 
critical value can be calculated |21( using the exact energy 
spectrum without linearization, yielding U c = 2.23|i|. 
The critical exponent of the order parameter s is [3 = 1 
which differs from the usual (thermal) MF value of 1/2. 
The difference can be understood from the expansion of 
the total energy for small s where the second order term 
is followed by a third order term (instead of ~ s 4 term 
which is usual in Landau theory): 
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with E(0) = Vf (1 

o-^2 . The unconventional third order term is induced by 
the linear behaviour of the DOS at the Fermi level. Near 
the critical point, 

E min (U) = E(0) -2N(U- U c f /(3W 2 ) 



(5) 

and the critical exponent of the ground state energy is 
a = —1. The control parameter in this quantum phase 
transition is U, as opposed to the reduced temperature in 
classical phase transitions. In an external magnetic field, 
we have E{s)/N = E(0)/N - as 2 + c\s\ 3 - hs close to the 
critical point, where h is the staggered field, coupled to 
the order parameter. At U = U c (a = 0), the extremal 
value of the energy is at \s\ = (h/(4W)) 1 / 2 so the critical 
exponent of the conjugate field is 5 = 2. Note, that for a 
normal Landau theory a — and 5 = 3. 

The bare susceptibility of the staggered magnetiza- 
tion at zero temperature is obtained as \ — 



dressed within RPA as 
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diverging at the critical point, and the critical expo- 
nent of the susceptibility is 7 = 1. The susceptibility 
can be determined as the inverse of the second deriva- 
tive of the total energy Q with respect to s evalu- 
ated at the minimum as well. This calculation provides 
X = (2a)- 1 = WiAU^iU-W)- 1 and also yields 7 = 1. 
The correlation length and the relaxation time can be 
written as £ = /wf/A and r = S/A so the corresponding 
exponents are v = 1 and z = 1. The critical exponents 
calculated above fulfill scaling relations. The renormal- 
ization group eigenvalues are yt = 1 and yh = 2 with 
the effective dimension d e ff = d + z = 3. The above 
exponents are summarized in Table |U 

Our results agree qualitatively with the RG analysis in 
the presence of long-range Coulomb interaction [22|, |23| . 
namely that around the Dirac point, weak interactions 
cannot trigger phase transitions. The critical proper- 
ties of the above MF theory for antiferromagnetism in 
graphene agrees with the BCS theory for graphene 24j, 
whenever comparison was possible. 
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TABLE I: Quantum MF exponents in graphene at T = 0. 

Having determined the MF quantum critical proper- 
ties of graphene, we observe that the exponents differ 
from those of a finite temperature thermal MF phase 
transition. To generalize these results, let us consider an 
electron system in d spatial dimension with energy spec- 
trum consisting of two bands, touching each other at the 
Fermi level with DOS 



P (e) = C\e\ r 



(7) 



where r > — 1 must hold so that the total number of 
electrons N does not diverge. We introduce a cut-off 
energy W so that N = CW r+1 /(r+l) holds, and the case 
r = 1 corresponds to graphene. Such a DOS follows from 
a non- interacting spectrum, £o(k) ~ ±fc z , connecting r 
with the dynamical critical exponent (z) as z(r + 1) = d 
By switching on an electron-electron interaction, an 
ordered phase can develop. Within the MF approxima- 
tion, an energy gap opens in the quasi-particle dispersion 
relation e(k) = ±^/e (k) 2 + A 2 . The control parameter 
of the transition is U and the order parameter can be 
denoted by s (describing arbitrary ordering with spin or 
charge or pairing character) , which is proportional to the 
energy gap A = Us. At zero temperature the total en- 
ergy can usually be written using the DOS and the energy 
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where NA 2 /U is the general energy cost of the ordered 
phase, 2-F\ is the Gauss hypcrgcomctric function. By 
expanding it with respect to A/W in order to determine 
the critical behaviour, we get 
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measured from the ground state energy in the normal 
(unordered) phase, E(0) = -NW{r + l)/(r + 2) and the 
coefficients are 
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where T is the gamma function. Higher order terms in 
the expansion are higher even powers of A/W (6, 8, ...), 
not influencing the critical behaviour. The ground state 
energy for r = 1 agrees with Eq. (U) for graphene apart 
from spin degeneracy. The expression of the energy, Eq. 
(|9]), differs from the usual Landau expansion because 
of the term A r+2 , producing different critical exponents 
when — 1 < r < 2. 

For r > 2, we get the usual MF critical exponents 
{a(U) changes sign and b > holds in this case), and the 
critical coupling U c = 2rW/(r + 1). The conventional 
Landau theory, with the temperature replaced by the in- 
teraction, operates. For — 1 < r < 2, the term A 4 can 
be neglected in © because close to the critical point it 
is much smaller than the A 2 and A r+2 terms. 

For < r < 2, the critical U is obtained from a(U c ) = 
which yields again U c = 2rW/(r + 1). Therefore, as 
long as r > 0, there is always a finite critical interaction, 
and no conventional weak-coupling solution exists. Due 
to this, near the critical point, a(U) ~ U~ 2 {U — U c ) 
so a(U) can be regarded as the deviation of the control 
parameter from its critical value, and its scaling can be 
used to determine quantum criticality. The total energy 
is minimal at 
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giving (3 = 1/r for a > and A = for a < 0. In the 
ordered phase the ground state total energy reads as 
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giving a = 1 — -. Upon introducing an appropriate ex- 
ternal field h which is conjugate to the order parameter, 
we obtain to linear order in small fields A = h/(2U\a\) 
in the normal phase (a < 0) and 8 A = h/(2Uar) in the 
ordered phase (a > 0). In both cases, 



X~ a 
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therefore 7=1. At the critical point a = 0, the total 
energy is extremal at 
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yielding 6 = r + 1. The exponents calculated above ful- 
fill scaling relations and yield the renormalization group 
eigenvalues y h = d e //^±|, y t = d eff7 ^ where d e ff = 
d + z = d(r + 2)/(r + 1) is the effective dimension of the 
system. 

This was the regime, often investigated in the pseudo- 
gap Kondo problem |19j. However, the — 1 < r < region 
is also possible physically 25| , and we focus on its criti- 
cal behaviour in the followings: c < and the coefficient 
a(U) is also negative for all U. The critical point is at 
U c = so we have to calculate with the original order 
parameter s, and not A. Near the critical point the total 
energy reaches its minimum at 
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leading to j3 = — 1 — - . The distinction between A and s 
is important when U c = 0, since A ~ [7 1 /' r ' J and its iden- 
tification as the order parameter would lead to incorrect 
scaling. The minimum energy scales as 
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determining a = 3+-. The linear response to an external 
field is Ss = —h/(2rU) and the susceptibility diverges 
with exponent 7=1, while a finite external field modifies 
the order parameter as s ~ h l+r ', giving S = 1/(1 + r). 

The scaling relations are fulfilled by these exponents 
and the renormalization group eigenvalues are y t = 
d e ff -^p; , Vh = d e ff -^ro ■ The correlation length exponent 
follows also as v = 1/yt- Note that the critical exponents 
calculated in the cases < r < 2 and — 1 < r < reach 
the same limit as r — >• (a — >• — oo, [3 —t oo, 7 — >• 1, 
5 — > 1), but are non-analytic at r — 0, as demonstrated 
below. The critical exponents for — 1 < r < are related 
to those for < r < 2 by replacing r in the latter by 
— r/(l + r), as seen in Fig. [T]and Table HU The energy 
gap scales with [/ 1 /' r L 

Surprisingly, the r = case requires separate consider- 
ations. The total energy acquires logarithmic corrections 
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close to the critical point, minimized by A = 2We u , 
so (3 — oo and ""^ - = —We~~u~, hence a = — oo. 
The correlation length scales as 1/A, implying v = oo. 
The exponential dependence of the order parameter and 
^ = 00 suggest that the conventional r = MF transition 
at T = as a function of the interaction and not temper- 
ature, resembles closely to a Kosterlitz-Thouless type of 
phase transition. Indeed, both the Kondo temperature 
in the Kondo problem and the gap in the sine-Gordon 
model scales similarly, but not identically to the above 
findings. The magnetic field dependence of the order pa- 
rameter reads as 



A, (A. 

W n \2W 



(19) 



defining (5 = 1 + 0, the "0" denoting the logarithmic cor- 
rections. Interestingly, the linear response susceptibility 
diverges as ~ 1/U 2 , giving 7 = 2. This non-analytic be- 
haviour is somehow expected, since all other exponents 
arc non-analytic at the origin, although their limiting 
values at r — > 0* agree. 




FIG. 1: The critical exponents as a function of the DOS expo- 
nent, r are plotted. For r > 2, the conventional MF exponents 
are recovered. 

In summary, we have studied the quantum critical 
properties of the MF transition in graphenc. An ordered 
phase requires a critical coupling to occur, and the crit- 
ical exponents listed in Table fl] differ from those in the 
Landau theory. The former suggests that weak interac- 
tion cannot trigger a phase transition in graphene, and 
can only renormalize some of its parameters. This ex- 
plains the lack of any experimentally verified phase tran- 
sition so far. Then, we have developed the MF theory 
of general bulk gapless phases, whose critical properties 
were found to vary strongly with the DOS exponent, sum- 
marized in Fig. Q] and Table |TTJ This extends the validity 
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TABLE II: Summary of the MF quantum critical exponents 
for general gapless phases. The exponents for — 1 < r < 
can be obtained from the < r < 2 exponents by the r — ¥ 
—r/(r+ 1) replacement. The r > 2 exponents are identical to 
the conventional thermal MF exponents. Mono- and bilayer 
grahene corresponds to r — 1 and with d — 2, respectively. 



of the MF universality class into quantum criticality, cov- 
ering a much broader region than that of the conventional 
Landau theory. 
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